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Global regularity for 2D Muskat equations with finite slope 

Peter Constantin, Francisco Gancedo, Roman Shvydkoy, and Vlad Vicol 

Abstract. We consider the 2D Muskat equation for the interface between two constant density fluids in an 
incompressible porous medium, with velocity given by Darcy’s law. We establish that as long as the slope of 
the interface between the two fluids remains bounded and uniformly continuous, the solution remains regular. 
The proofs exploit the nonlocal nonlinear parabolic nature of the equations through a series of nonlinear lower 
bounds for nonlocal operators. These are used to deduce that as long as the slope of the interface remains uni¬ 
formly bounded, the curvature remains bounded. The nonlinear bounds then allow us to obtain local existence 
for arbitrarily large initial data in the class 1 < p < oo. We provide furthermore a global regularity 

result for small initial data: if the initial slope of the interface is sufficiently small, there exists a unique solution 
for all time. July 7, 2015 
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1. Introduction 


In this paper we study solutions of the one dimensional nonlinear nonlocal equation 

5af{x,t)a 


d, 


= f 

JR 


{5af{x,t)Y + 


da, 


where we have denoted by 


daf{x, t) = f{x, t) - f{x - a, t) 
the a-finite difference of / at fixed time t > 0, and by 

f{x,t) = dxf{x,t) 

the partial derivative of / with respect to the position a: G M. 
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1.1. Physical origins of the problem. Equation (1.1) is derived from a two dimensional problem, 
where y = f{x, t) describes a graph (in coordinates whereby gravity acts in the y direction) which represents 
the interface separating two immiscible fluids of different constant densities. The flow is incompressible, 
and the fluids permeate a porous medium. The 2D Darcy’s law [Dar56] 

^u{x,y,t) = -Vp{x,y,t) - {0, gp{x,y,t)) 

is used to model the balance of gravity, internal pressure forces and velocity (replacing acceleration). Here u 
and p are the velocity and pressure of the fluid, which are functions of time f > 0 and position (x, y) G M^. 
The physical constants of the model are viscosity p, permeability k, and gravity g. For simplicity we fix 
these physical parameters to equal 1. We pay special attention to the density function p, which is given 
by a step function representing the jump in the fluids densities. Taking the size of this jump equal to 27r 
(denser fluid below), we arrive at (1.1) (see [CG07] for a detailed derivation). This is the classical Muskat 
problem [Mus34], which has been broadly studied [Bea72] due to its physical relevance, wide applicability, 
and mathematical analogy to a completely different phenomenon: the interface dynamics of fluids inside 
Hele-Shaw cells [ST58]. 

The Muskat problem has been widely studied, its multiple features being taken under consideration. 
It can incorporate capillary forces that deal with surface tension effects [CP93, C93, KM15], density 
jumps [CG07], viscosity [SCH04], and permeability discontinuities [BCGB14]. The dynamics is different 
for the one fluid case (with a dry region) [CCFG13a], two fluids [SCH04], and multi-phase flows [CGIO]. 
Moreover, one may consider boundary effects [GB14] and study the interface in the three dimensional 
case [A07, CCG13b]. 

1.2. Summary of known results. The basic mathematical questions regarding (1.1) are the existence 
and regularity of solutions. The Muskat problem can be ill-posed, due to existence of unstable configura¬ 
tions [Ott99, Szel2]. This feature can be understood in the contour dynamics description via the Rayleigh- 
Taylor condition, which involves the geometry of the density jump or the dynamics of fluids with different 
viscosities. Essentially, the unstable situation occurs when a denser fluid is above a less dense fluid, or when 
a less viscous fluid pulls a more viscous one [A07]. Surface tension effects regularizes the equation [ES97, 
A14], so that the system is well-posed without satisfying the Rayleigh-Taylor condition, but there are still 
instabilities in this case [GHS07, EMU]. In the stable cases it is possible to find global-in-time solutions 
for small initial data [SCH04, CG07, CGBS14] due to the parabolicity of the equation. This smallness 
can be measured in different regularity classes and sizes [CCGS13, CCG+13a, GB14, BSW14, CGBS14]. 
Eow regular global existence results are obtained for Eipschitz weak solutions in the case of small slopes in 
[CCGS13, CCG+13a]. 

A very interesting phenomenon for the Muskat problem is the development of finite time singulari¬ 
ties starting from stable initial data. One possible singularity formation happens by the breakdown of the 
Rayleigh-Taylor condition in finite time [CCE+12]. In this case the interface cannot be parameterized by a 
graph of a function any longer, i.e. there exists a time > 0 when solutions of equation (1.1) satisfy 

lim ||/'(f)||Loc = oo. 

t—¥t\ 

After that time the free boundary evolves in a way that produces a region with denser fluid above less 
dense fluid. For short time the interface is still regular (although it is not a graph) due to the fact that 
the parabolicity of the system gives instant analyticity. This sequence of events is termed “wave-turning” 
because it is a blow up of the graph-parameterization. Wave-turning has been proven to arise in more general 
cases with different geometries [BCGB14, GSGB14]. The regularity of the interface is lost at some time 
t 2 > ti after wave-turning, and the interface ceases to belong to C^. Thus, starting from a stable regime, 
the Muskat solution enters an unstable situation, and then the regularity breakdown occurs [CCFG13a]. 
This phenomenon is far from trivial, as some unstable solutions can become stable, and again reach unstable 
regime [CGSZ]. 
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A different type of singularity which may develop in solutions of the Muskat problem is given by the 
occurrence of a finite time self-intersection of the free boundary. If the collapse is at a point, while the 
free boundary remains regular, then the singularity is termed a “splash”. If the collapse is along a curve, 
while the free boundary remains regular, then it is termed a “splat” (or “squirt singularity”). Both kinds of 
collisions have been ruled out in the case of a stable density jump in [CGIO, GS13]. In the one fluid case, 
there exists a finite time splash blow up [CCFG13b] in the stable regime, but the splat singularity is not 
possible [CPC14]. 

1.3. Main results. The main purpose of this work is to develop conditional regularity results (blow 
up criteria) for the Muskat problem involving the boundedness of the first derivative of the interface. Our 
motivation is both to complement and to aide the theory of singularities for the Muskat equation by obtaining 
sharp regularity results in terms of purely geometric quantities such as the slope. 

The bounds of nonlinear terms used to obtain the conditional regularity results also allow us to obtain 
several existence results in low regularity regimes. Specifically, we prove local-in-time well-posedness 
and global well-posedness for Muskat solutions in classes of functions of bounded slope and (with 
p € (1, oo]) integrability of curvature. 

Recently, results involving solely the second derivative of / appeared in the work [CGBS14], where 
the authors explore the PDE structure of the original Darcy’s law in the bulk rather than the corresponding 
integral equation for contour dynamics (1.1). They develop an well-posedness theory for the equation 
in the case of density and viscosity jumps. Both the local and global existence results are obtained under 
smallness assumptions for the interface: iT3/2-i-e local results and for the global results. These 

smallness conditions imply /' E for some /? > 0, requiring thus smallness of the Holder continuity of 
the slope. 

The optimal well-posedness theory (local, conditional, or global) for the Muskat equation should (con- 
jecturally) involve only assumptions of uniform boundedness of the slope /'. To the best of our knowledge, 
this remains an outstanding open problem. The difficulty in reaching a regularity theory for / may 

be seen in at least two ways: 

• When considering the evolution of /', linearized around the steady flat solution, the resulting equation has 
L°° as a scaling invariant norm. Therefore, for large slopes, obtaining higher regularity of the solution 
(required in order to obtain uniqueness) involves solving a large data critical nonlinear nonlocal problem. 
Of course, the additional complication is that, due to the fundamentally nonlinear nature of the equation, 
in the large slope regime the linearization around the flat solution rapidly ceases to be useful, and new 
severe difficulties arise. 

• The velocity v transporting /, /', and /" (cf. (2.2), (2.5), and (2.6) below) is obtained by computing a 

Calderon commutator applied to the identity (cf. (2.3) below). However, when /' is merely bounded, 
while the Calderon commutator maps —)• for p E (1, oo), the endpoint L°° inequality fails [MS13]. 

Moreover, since we are in ID the x-derivative of v enters the estimates (we cannot use incompressibility), 
and this term (cf. (2.11)) yields a Calderon commutator applied to 6af'{x,t)/a. Thus, to bound this 
requires estimates on the maximal curvature, not slope. This is the reason why up to date the existing 
continuation criteria [CG07, CGBOI14] required that the solution lies in for some e > 0. 

In this paper we move closer to the conjectured critical / E W^’°° well-posedness framework. The 
main new idea is that the maximal curvature can be controlled by either assuming uniform continuity of 
the slope /' (cf. part (i) of Theorem 1.2), or by assuming that the initial slope ||/ol|oo is sufficiently small 
(cf. part (ii) of Theorem 1.3). We moreover show that under either of these conditions one can also control 
||/"||iP for all p E (1, oo). In the regime p E (1, 2) not even the local existence of solutions was previously 
known. 

The important novelty of our large data result is that the maximal slope (and the maximal initial cur¬ 
vature) are allowed to be arbitrarily large, as long as the slope obeys any uniform modulus of continuity. 
Our uniform assumption on /' is a local smallness of variation assumption instead of a global smallness 
assumption of /' in L°°. The modulus of continuity of /' is not required to vanish at the origin at any 
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particular rate, and it can be even weaker than a Dini modulus (and thus the Calderon commutator in v and 
the nonlinear terms would not necessarily yield bounded functions). 

Regarding our small data result, we draw the attention to the fact that the smallness is assumed only on 
the initial slope, not on the W‘^’P norm. 

The main tools used in our analysis are various nonlinear maximum principles for the evolution of /", in 
the spirit of the ones previously developed in [CV12, CTV15] for the critical SQG equation. The robustness 
of the pointwise and integral lower bounds available for the nonlocal operator Cf defined in (2.8) below 
allows us to treat all values of p E (1, oo] and enables us to analyze the long time behavior of the curvature, 
cf. (1.5). The endpoint case which scales as remains however open. 

Next we state our results more precisely. First we establish a low regularity local-existence result for 
(1.1) with initial datum that has integrable, respectively bounded, second derivatives. 

Theorem 1.1 (Local-existence in Assume the initial datum has finite energy and finite slope, 

that is /o E n Let p E (l,co], and additionally assume that /o E Then there exists a 

time T = r(||/o||w 2 ,pntycoo) > 0 and a unique solution f E L°°(0, T; fl (7(0, T; n 1F^’°°) of the 
initial value problem (1.1) with initial datum f{x, 0) = foix). 


The main result of this paper is: 

Theorem 1.2 (Regularity criterion in terms of continuity of slope). Consider /o E for k > 3 

and assume that f is bounded on [0, T], i.e. that 

sup ||/'(f)||Loo < .B < oo (1.2) 

ie[o,T] 

for some B > 0. Furthermore, assume that f is uniformly continuous onM x [0,T]. That is, assume there 
exits a continuous function p: [0, oo) [0,oo), that is non-decreasing, bounded, with p{0) = 0, such that 
f obeys the modulus of continuity p, i.e. that 

\Saf{x,t)\ < p{\a\) (1.3) 

for any x, a E M and t E [0, T]. Then the following conclusions hold: 

(i) SUPjg[0,T] ll/"(t)l|L°° <C{\\fl^\\L^,B,p). 

(ii) The solution stays regular on [0, T] and f E (7([0, T]; H^) n L^(0, T ; ). 


We note that in Theorem 1.2 we do not require the modulus of continuity p to vanish at 0 at any specific 
rafe. For example, if can be guaranfeed by a condifion /' E L°°(0,T;(7^) for /3 > 0. The remarkable 
fealure of conclusion (i), however, is fhaf fhe confrol on /" is furnished by fhe inifial dafa only. If however 
fhe inifial slope /' is small enough, fhen we can show fhaf fhe condifion f" E can be propagafed wifhouf 
any furfher assumpfions on fhe confinuify of fhe slope. And fhaf way we can obfain a global existence in fhe 
corresponding class. 

Theorem 1.3 (Small slope implies global existence). Consider initial datum /o E Lf, with maximal 
slope that obeys 

\\fo\\L^<B<^, (1.4) 

for a sufficiently large universal constant ( 7 * > 1. Let p E (l,co]. If we additionally have /q E L'P, then 
the local in time solution obtained in Theorem 1.1 is in fact global, and we have 

||/"(f)||LP < m\Lr 


for all f > 0. Furthermore, we have that 

(i) for p = oo the curvature asymptotically vanishes as 


II/," 


0 \\L° 


1 + 


II/o 


1005 




(1.5) 
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for all t > 0, 

(ii) for p E [2, oo) the solution obeys the energy inequality 


wnmi, + + 3„()b(i\b2) I‘ iirwiiK-;.* £ iirmiiw d-e) 


\LP '1 + ^2 
for all t > 0, 

(Hi) for p E (1, 2) the solution obeys 


1 


ii/"(^)ii 


r/s < 

LP+1^^ \ 


CB{l + B^) 7o 

for all t > 0, for some sufficiently large constant C > 0 that may depend on p. 


p 

LP 


(1.7) 


We note that the assumption (1.4) together with the maximum principle for /' established in [CG09, 
Section 5] show that (1.4) holds at all times f > 0. Now we see that bounds established in (i), (ii), (iii) also 
prove maximum principles at the level of the curvature. We complement existence result with the following 
uniqueness statement which in part shows that solutions of Theorems 1.1, 1.2, 1.3 are unique. 


Theorem 1.4 (Uniqueness). Let /o E L°°, and assume that f is a solution to (1.1) on time interval 
[0,T] with the following properties: sup^gjQ'j’] ||/^(f)||L°° < oo; f{x,t) 0 as |x| ^ oo, for all t < T; 
dtf{x, t) exists for all (x, t) and ||5i/||loo(kx[0,t]) < oo; /' is uniformly continuous on M x [0, T], i.e. (1.3) 
holds for some modulus of continuity p. Then f is the unique solution with fo{x) = f(x,0) satisfying all 
the listed properties. 


We note that the only essential assumption in the list of Theorem 1.4 is the same uniform continuity 
from Theorem 1.2. The remaining assumptions are present to justify the application of the Rademacher 
Theorem, see the Appendix, and we believe these additional assumptions may be avoided. For instance, a 
sufficient condition which ensures that dtf{x, t) exists for all (x, t) is that /" E Lf^Wx'^ for some p > 1, 
or that the modulus p obeys the Dini condition. 

We also note that in the course of proving Theorems 1.1 and 1.3 we follow the standard strategy: we 
obtain all necessary a priori estimates and construct solutions by passing to the limit in the regularized 
system as elaborated for instance in [CCG+lSa]. To avoid redundancy, the proofs in this paper only consist 
of these a priori estimates. 


2. Preliminaries 


For the remainder of the manuscript we shall use the following notation for finite-difference quotients: 

^af{x) 


A«/(x) = 


a 


for X E M and a E M \ {0}. After integration by parts (1.1) becomes 

f{x,t)a - 6af{x,t) 

((5a/(x,t))2 + a2 

Saf{x,t) 


dtf{x,t) = p.v. [ 

Jr 


da 


or equivalently, 

dtf{x,t) + v{x,t)dxf{x,t) + P.V. [ 

Jr 

where we have defined fhe velocify field 


(dafixp))"^ + 


da = 0 


v{x,t) = —P.V. 

Jr 


a 


rda = -P.V. 


1 


da 


(A«/(x,f))2 + l a 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


((5Q,/(x,f))2 + 

In particular, from (2.2) if is immediafe fhaf fhe maximum principle for fhe global bound on / holds (see 
[CG09]). Moreover if is known fhaf fhe maximum principle also holds for fhe energy (see [CCG’'"13a]), i.e. 


wfmip < wmilp 


( 2 . 4 ) 
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holds for p = 2 and p = oo. 

2.1. Equation for the first derivative. The equation obeyed by f'{x, t) is 

6 af{x,t) 


dt{f)+v{x,t)dM')+Py- f 

Jm 

= 2 / 


{ 6 afix,t)f + 
i 6 af{x,t) - af{x,t)){5afix,t)){6af'ix,t)) 


da. 


(((5a/(x,t))2 + 0(2)2 

2.2. Equation for the second derivative. The equation obeyed by the second derivative is 


(2.5) 


dtin+v{x)dM'')+py- [ 

Jr 




rda 


= 4P.V. 

Jr 

+ 2P.V. 
+ 2P.V. 

- 8 P.V. 


i^afixW + a'^ 

{Saf{x) - af"{x)) 6 af{x, t) 6 af{x, t) 


da 


da 


ii^afix))^ + a2)2 
i^afix) - af'ix)){6af{x)f 

{{^afixW 

i^afix) - af'{x)){Saf{x)){Saf'{x)) 

((5„/(x))2 + a 2)2 

{^af{x) - af(x)){5afix))^{6afix))^ 


{{^afixW + a'^f 


da 

da 


=: T1+T2+T3 + T4. 

We now pointwise multiply (2.6) by f"{x, t) and obtain 

{dt + vd, +Cf) |/"(x, t)|2 + Vf[f]ix,t) = 2fix,t) (Ti + Ta + Ts + T 4 ) 
where we have denoted 

da 9 {x) 


and 


Cf[ 9 ]{x) = P.V. [ 
Jr 

Vf[g]ix) = P.V. [ 
Jr 


((5«/(x))2 + a2 
iSa9ix)f 


da 


da 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 


(<5a/(x))2 + 0(2 

for any smooth function g, and Ti,..., T 4 are as given by (2.6). It will be useful to also consider (2.7) where 
the transport term is written in divergence form, namely, 

(dt +£f) Iff + dMn^) + Pflf"] = V" {Ti + Ta + Ta + T 4 ) + \ff Tg 

where we have denoted 

6afix,t) Aafix,t) 


T^{x,t) = dxv{x,t) = 2P.V. [ 

Jr 


a2 (l + (A„/(x,f))2)2 
Note that when f = c, where c is a constant, the above operators become 

6 ag{x 


da. 


( 2 . 10 ) 

( 2 . 11 ) 


C,[g]{x) = P.V. [ 
Jr 

V,[g]{x) = P.V. 


a^ 


-da =: Ag{x) 


{Sa9{x))‘^ j ^ 

-^- da =: T>[g\{x). 




We further note that 


[ P[ 9 ]{x)dx = 
Jr 


|2 

1 + 1/2 


( 2 . 12 ) 


(2.13) 
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and 


'Df[g]{x) > 


(2.14) 


for all X E M. 


2.3. Evolution of the pth power of the second derivative, for p >2. When p E [2, oo), we consider 
the function 

cpp{r) = |r|P/^ 

and we multiply (2.7) by ip' {\f"{x, f)P). The convexity of (pp (which holds if p > 2) ensures that 


p^p{\nx,t)\^)Cf[\nx,tr] > Cf[ppi\nx,t)\^)] 


and we thus obtain 


{dt + vd, +£f) |/"(x, t)\P + P\n^^^}^Vf[r]ix,t) < p|/"(x, t)r^ iTi +T2 + Ts + T^\ . (2.15) 
Moreover, we may again write this in divergence form as 


{dt+Cf)\np+dMf"\n + 

where is as defined in (2.11). 


%[/"] < p\fT-^ iTi +T2 + T3 + T^\ + l/'r l^sl ( 2 . 16 ) 


2.4. Evolution of the pth power of the second derivative, for p E (1,2). When p E (1,2), the above 
trick does not work, since the function pp is in this case concave. Instead, we consider the function 

^p(r) = \r\P 

which is convex, and obeys 


jr-s? 


'ippis) - il^pir) + 'iljp{r){r - s) = |s|p - \r\P + p'-^{r - s) > p(p - 1) \ 


4 ’ 


2|r| > |s| 
Isl > 2|r| 


(2.17) 


for all r, s E M. Thus, for p E (1, 2) we are lead to consider the nonlocal operators 




and 


[9]{x) = [ 

Jr 


V^[ 9 ]{x)= [ 
Jr 


\g{x) - g{x - a)\P 
o? + (<5a/(x))2 

\g{x) - g{x - a)\P 


da 


da. 


Note that for p = 2, we have V 2 = V, and that for g E n the above integrals are well-defined 
(even without a principal value), since p > 1. Using (2.17) we may show that upon multiplying (2.6) by 


'ipp{f”ix,t)) =p|/"(x,f)|P//"(x,f), that 


{dt + v{x,t)dx + Cf)\f"{x,t)\P + 


P(P-I) . / 2 T>/[/"](a^>i) 


• mm 


< p|/"(x,f)|P"^|Ti(x,f) + . . . +r 4 (x,f)|, 


(2.18) 


or in divergence form. 


{d, + £f)\fr + dMf"n + min 


<p\f"\P-^\T, + ... + T^\ + \f"\P\n\, 


(2.19) 


in analogy with (2.15) and (2.16). 
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2.5. Bounds for Taylor expansions. The finite difference quotient may be bounded directly as 

|A„/(x)| < B (2.20) 

for any a, x. We note that we may expand 

7^l[/'](x, a) = Aaf{x) - f'{x) = - f if{z) - f{x))dz. 

^ Jx—a 

The bound 

|7^l[/'](x,a)| <25 

is immediate due to ||/'||£,oo < B, but if additionally /' has a modulus of continuity p, then we have the 
improved bound 

\Aaf{x) - f{x)\ < p{\a\). (2.21) 

Without loss of generality we assume p is not linear near the origin, since then the conclusion and the 
assumption of the theorem are identical. 

It will be convenient to denote 


7^l[/"](x,a) := Aaf'{x) - f''{x) 

a I 


inz) - nx))dz 


( 2 . 22 ) 


and 


a 


T^ 2 [f']{x,a) := Aaf{x) - f{x) + ^/"(x) 


= 7^l[/'](x,a) + |/"(x) 


- r [\f"H - f"ix))dwdz. 

® Jx—a J X 


(2.23) 

/ x—aJX 

as the first order and the second order expansions of f''{x — a) around f"{x). For these terms we have 
pointwise in x estimates and a in terms of the dissipation present on the right side of (2.7). First, we have 
that 

|/"(z) -/"(x)| 


|7^l[/"](x,a)|<- r 

^ J X- 


Z — X 


\z — x\dz 


a 


\z — xfdx 


\ 1/2 


(2.24) 


for any a > 0, for some universal constant C > 0. The bound for all a / 0 trivially holds upon replacing 
with Similarly, it follows that 


X f'Z I XU 


iw2i/"i(x,o)i<i r / 

^ J X — OL J X 

If (f 

^ J X—a. \j X 


\f"[w)-f"{x)\ 


— x| 


w — x\dwdz 
dw 


< C 


\w — X 


\ 1/2 

Im — x\‘^dw j dz 


(P|/"](x))l /2 


a 


f \z-xf^^d2 

J X—a 




(2.25) 
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for any a > 0, for some universal constant (7 > 0. The bound for all a / 0 trivially holds upon replacing 
0 -3/2 |q;| 3/2 Lastly, we note that for p G (1,2) in a similar way to (2.24) and (2.25) we may bound 

|7^l[/"](x,a)| < CaVP(pP[/"](^))i/p (2.26) 

|7^2[/"](x,a)| < (2.27) 

with a constant C = C(p) > 0. 


3. Nonlinear lower bounds 


In this section we use (1.2) in order to obtain lower bounds for Py[/"](x) at any x G M. These lower 
bounds follow in a similar way to the bounds previously established in [CV12, CTV15] for A. The main 
results of this section are. 

Lemma 3.1. Assume that f is Lipschitz continuous with Lipschitz bound B, i.e. that (1.2) holds. Then 
we have that 

»/l/"l(^) > n.o-/. 0-1) 


(3.2) 


245(1 + ^ 2 )' 

holds pointwise for x G K. Moreover, for p G (1, 2), we have that 

holds for all x G M. 

The above lower bound however will only suffice to prove a small B result. For large values of B we 
must obtain a better lower bound. If we additionally use a modulus of continuity obeyed by /', we obtain: 

Lemma 3.2. In addition to the assumption in Lemma 3.1, assume that f obeys a modulus of continuity 
p. Then there exists a continuous function Lb ■ [0, (X)) [0, (x) that implicitly also depends on p, such 

that 


Vf[f"]{x)>LB{\f"{x)\) 
holds pointwise for x G M, and we have that 

r ^siv) 

lim -^— = cx 

y^oo y'^ 

at a rate that depends on how fast lim^^g+ p[r) = 0 . 


(3.3) 

(3.4) 


Next we provide two lower bounds whose proofs follow by similar arguments to those in Lemma 3.1 
(see also [CV12, CTV15]). 

Lemma 3.3. Assume that f is Lipschitz continuous with Lipschitz bound B, and let p G (1, cx). Then 
the following lower bound 

holds for all x G M. Moreover, for p G (1, 2) we have that 

holds. 


(3.5) 


(3.6) 


Lemma 3.4. Assume that f is a regular function with Lipschitz bound B. Then the following lower 
bound holds for any x G M. 

1 .o 


2A\\f"\\L^{l + B‘^) 


'DAHix) > 




(3.7) 
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Proof of Lemma 3.1. For this purpose, let r > 0 to be chosen later, and let x be an even cutoff 
function, with x = 1 on [1, oo), x = 0 on [0, 1/2], and x^ = 2 on (1/2,1). We then have 

f \f"{x)\^-2f"{x)f"{x - a) /O' 




> 


1+52 
1 


1+52 

Integration by parts yields 




(1+^2) 


1 + .B2 

/ da 

2|/"(x)|2 


a^ 


X - da 


2|/"(x)| [ 

Jr 


da fix — a) /a 
-X 




r 

da 


-4|/"(x)| 


|/'(x-a)| 


-4|/"(x)|- 


da 

\a\>r/2 

f \fix-a)\ 

Jr/2<\a\<r a 


da 


> 


1 /2|/"(x)|2 24i?|r(x)| 


( 1 +^ 2 ) 


Letting 


2AB 


r = 


\f'ix)\ 

in the above estimate, we arrive at the bound (3.1). 

Let p E (1, 2). In order to prove (3.2) we appeal to the inequality 

|r — > jrj^ (1 — p- 


which holds for any r, s E M and all p > 1. It follows that 

a 2 tTIp L 


> 


1 + ^2 

1 


\f'ix) - fix - a)\P (a 
-X 




(1 + Bf 


\fix)\P p|/"(x)|P-i 

a2 (1 + ,B2) 


da 


da fix — a) fa 
0 X 
a^ 


da 


> 


1 \fix)\P 24 B\fix)\P 


-1 


1 + .B2 


Upon choosing 


1 + .B2 


r = 


485 


the proof of the Lemma is completed. 

Proof of Lemma 3.2. Similarly to the above proof, we bound Vj from below as 

1 / 2 |/"( x )|2 


^finix) > 


> 


> 


> 


(1 + 52 )' 

1 ( 

V r 

.2ir(x)|2 

(1 + 52 )' 

V r 


<2|/"(x)|2 

(1 + 52 )' 

V r 


.2ir(x)|2 

(1 + 52 )’ 

V r 


-2ir(x)i [ 

Jr 

-4|/"(x)| ' 

-4|/"{x)|i f 

r Jr 

- 8|/"(x)| [ 

Jr^ 

-12|/"(x)| / 

J r 


da if ix - a) - f ix)) (a\ 

-X [-] da 


a^ 


\^afix)\ 


da 


|Q|>r/2 1^1 


|d«/'(x)| 


r/2<\a\<r ^ 

pia 


da 


r/2 a' 


3 

pia 


da — 

da 


1 r 

8|/"(x)|- / 

P Jr/2 


pia) 


da 


a^ 




(3.8) 


(3.9) 


□ 
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In this case we choose r = r(|/"(x)|) as the smallest r > 0 which solves 


\nx)\ 

6 


= r 


^oo 

Jr/2 


p{a) 


da. 


a-^ 


(3.10) 


The existence of such an r is guaranteed by the intermediate value theorem, and by computing the limits at 
r —0+ and r —oo of the right side of (3.10). Indeed, these limits are 


lim r 
r—^0 


I 


p{a) 


since p is not Lipschitz, and 


Note that with this choice we have 


lim r 

r^OQ 


r/2 

p{a) 


da = +00 


poo 

Jr/2 


da = 0. 




lim r(y) = 0. 

y^oo 

With the choice (3.10), the lower bound obtained on 'Df[f''] becomes 

P/[/"](x)>Ls(|/"(x)|), 

where the function Lsiy) is defined on (0, oo) implicitly by 


Lsiy) = 


with 


Jy) f 

J r\ 


p{oi) _y 

Q /? ‘ 

r{y)/2 ^ ^ 


(3.11) 

(3.12) 

(3.13) 


(1 + 52)r(y)’ 

A short computation shows that since r{y) is Lipschitz continuous on (0, oo), which implies that is also 
Lipschitz continuous on this interval. The important aspect to notice is that since /9(0) = 0, we have 


r Lsiy) 

lim -=— = oo, 

y^oo 

which means that the lower bound (3.12) is sharper than (3.1), when \ f"{x)\ S> 1. 
In order to prove (3.14), note that 

LB{y) 1 


(3.14) 


?/3 (1 + B^)yr{y) 


and 


yrjy ) 
6 


= r{yf j 

J r\ 


p{a) 


da. 


r{y)/2 




Therefore, in view of (3.11), the limit in (3.14) indeed diverges if we establish that 

re^Ja = 0. 

r^.0+ Jr/2 a 

Clearly, it is sufficient to verify that 

Ito ,2 /■'^ = 0. 

r^0+ Jr 

Indeed, since /? is a modulus of continuity of /' and we have that ||/'||ioo < B, we obtain 


(3.15) 


poo 

7 

JyJF 


p{a) 


da < 2Br^ 


a'^ 


roc 

J^ 


1 


da < Br ^ 0 


a 


3 — 


On the other hand, it is easy to see from the monotonicity of p that 


rVr 

1. 


p{a) 


a-^ 


da < p{y/r)r‘^ 


da piVr) 


as 


0 as 


r —)• 0. 


a-^ 


since p{0) = 0. 


□ 
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Proof of Lemma 3.3. In order to prove (3.5) we proceed as in (3.8), but we do not integrate by parts 
in the second term. This yields 




|/"(x)|^ 


> 


+ 52 r 

2 ir(x)p 

1 + 52 j- 


:\nx 


1 


{p-l)/p 


LP 


1 + 52 

8 \nm\n\LP 

1 + ^(p+i)/p 


/|a|>r/2 |a|2P/(P-l) 


for all p > 1. The desired inequality follows upon choosing 

nf'TLP 

\f"{x)\P 


Similarly, for p G (1, 2) in order to prove (3.6), we consider (3.9), but we do not integrate by parts in 
the second term. We obtain 

^. 1 _ 4 \nx)rm"h. 

^ + r 1 + ^2 r(P+i)/P 

Choosing exactly as above concludes the proof of the Lemma. □ 



Proof of Lemma 3.4. We consider the inequality analogous to (3.8) with /" replaced by /'". After 
integrating by parts once, the same argument used to prove Lemma 3.1 yields (3.7). □ 


4. Bounds for the nonlinear terms 

In this section we give pointwise in x bounds for the nonlinear terms Ti{x), with z G {1,..., 5} appear¬ 
ing on the right sides of (2.7) and (2.10). We first fix a small constant 

£ G (0,1] 

to be chosen later. The bounds we obtain depend on this e, and e will be chosen differently in the proofs of 
Theorem 1.1, Theorem 1.2, and Theorem 1.3 respectively. The main result of this section is: 

Lemma 4.1. Let B > 0 be such that (1.2) holds and fix e G (0,1]. There exists a positive universal 
constant C > 0 such that the bounds 

|ri(x)| + \T2{x)\ + |r3(x)| + |r4(x)| < ^|/"(x)|2 + (4 .i) 

|r5(x)| < CB\Anx)\ + ^|/"(x)| +C£2i?2^^^ (4.2) 

hold for all X GM. where f"{x) 0. 

In addition, we need a pointwise estimate for the nonlinear terms in terms of the dissipative term T>P[f"], 
whenp G (1, 2). 

Lemma 4.2. Let B > Qbe such that {1.2) holds, letp G (1, 2), and fix e G (0, 1]. There exists a positive 
universal constant C > 0 such that the bounds 

\TVx)\ + \T 2 {x)\ + |r3(x)| + |r4(j,)| < -5 ^|/''(x)| 2 + CeB"2!lG|W. ( 43 ) 

|r3(x)| <CB|A/'(x)l + ™|/"(x)|+CeB^2d^ (4.4) 


hold for all X GM. where f"{x) 7 ^ 0. 
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Proof of Lemma 4.1 . Throughout this proof, we fix a cutoff radius (for a) 

sB 

rj = rj[x) = 


\f {x)\ 

Note that at the points x where f''{x) = 0, there is no estimate to be done for Ti{x) terms as they are 
multiplied with f'{x) in (2.7) and (2.10). 

Estimate for the Ti term. We decompose T\ into an inner piece and an outer piece according to 

(Safix) - af"{x)) Safix) 6af'ix) 


Ti{x) = 4P.V. 
= 4 P.V. 

+ 4 


7^l[/"](x,a) A«/(x) (/"(x)+7^l[/"](x,a)) 


J\a\<rj ((A„/(x))2 + 1) a 

f T^i[f'']{x,a) Aaf{x) 6af'{x) 

’\a\>ri 


da 


-da 


((A„/(x))2 + l)^a2 

=: Ti^inix) +Ti^out{x). 

Using the pointwise in x and a bounds (2.20) and (2.24), and the definition of rj in (4.5), we obtain 


|ri,,„(x)| <Ci7|/"(x)|(P[/"](x))V2 


da 


'H<v 


\a 


11/2 


+ CBV[f"]{x) 


da 


'\a\<r] 


< Ce^/^B^/^\f"{xtl\V[f"]{x)f/^ + 


< CB\f{x)\^ + C£B 


2 nn{^) 

\f''{x)\ 


(4.6) 


and 


|ri,<n.,(i)l < CB\v[f%x)Yi-^ [ ^ 

J\a\>ri |Q:| 

< I /"(x) I {V[f'] (x)) 1/2 

for some universal constant C > 0. Combining (4.6) and (4.7) we arrive at 

2 |/"(x)ri(x)| < Ce-^B\f"{xt + CeB‘^V[f]{x) 

for some universal C > 0. This bound is consistent with (4.1). 

Estimate for the T 2 term. We decompose T 2 as 

(Aa/(x) - f{x)){6af'{x)f 


(4.7) 

(4.8) 


T2 = 2 P.V. 
= 2 P.V. 

= 2P. 


((A„/(x))2 + 1)2^3 

{-af"{x)/2 + 7^2[/"](x, a))((5a/'(x))2 


da 


'.V. I 

J\c 


((A„/(X))2 + 1)2«3 

7^2[/"](x,a)(5„/'(x))(r(x)+7^l[/"](x,a)) 


H<v 


- fix) P.V 


((A„/(x))2 + 1)2^2 

(/"(x)+7^l[r](x,a))2 


da 


da 


+ 2 


J\a\<ri ((Aq,/(x))2 + 1) 
(-a/"(x)/2 + 7^2[/"](a;, a))(5a/'(x))2 


J\a\>n ((Aa/(x))2 + l)2a3 

— • 2^2,1,in “t“ ^2^2,in “1“ 2^2,ont* 


da 
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By appealing to (2.25) and (2.24) we may estimate the inner terms as 


=))■'■" fi 


|r2,i,i„| < CB(vinx)rr‘ I |/"(i)| f ^ 


+ j da 

J\a\<r) , 




< CB\f{x)Y + CeB 


.2 nn{^) 


(4.9) 


and 


\T2,2,in\<C\f''{x)\i\f''{x)\^ [ da + \f"{x)\{V[f]{x)f/‘^ [ \a\^/‘^da + V[f"]{x) [ \a\da 

\ ^|q ;|<?7 J\oi\<rf 


\a\<r} 


< CB\f"(x)\^ + 

while the outer terms may be bounded as 


'|o|<?7 


(4.10) 


\T2„A < CB^ (\r{x)\ f ^ ^ 

Y J\a\>r] 1^1 J\a\>r} |^| ^ 

< Ce-^B\f{x)f + Ce-^/^B^/^\f{x)\^/^{V[f]{x))^/‘^ 


< Ce-^B\f{x)\'^ + CeB 


2nnM 


(4.11) 


for some universal C > 0. Combining (4.9), (4.10), (4.11), using the identity (2.13) and the Cauchy- 
Schwartz inequality we arrive at 


2 \fix)T2ix)\ < Cs-^B\r{x)\^ + CsB^Vinix) 

for some universal C > 0. This bound is consistent with (4.1). 

Estimate for the T 3 term. We bound T 3 as 


1^31=2 


py f - f{x)) Aa/(x) 6 af"{x) ^^ 


a 


((A„/(x))2 + 1)2 a 


<CB{V[f"]{x))^/^ ( [ (^«/(^) ‘ 




<CB{V\f'\{x)fl‘^ 


(a/"(x)/2 + 7^2[/"](x,a]) 


2\1/2 


'|o|<»7 




+ Ci?(P[/"](x))^/2 


(7^l[/'](x,a)) 


2\V2 


'|o|>'7 




—: T^,in + T3 


out' 


(4.12) 
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Using (2.25) we may further estimate 


T 3 ,in < CB{V[f"]{x))^/^ 
< CB{V[f"]{x)y/^ 


'H<v 


ianx)/2 + n2[f"]ix,a]))^Y" 

a2 


X f l/"(a;)l^ / da + \f''{x)\{V[f']{x))^^‘^ \a\^/‘^da + V[f”]{x) 


1/2 




'\a\<r] 


'\a\<-n 


\a\da 


< 


< C’SeV25l/2|j"(^)|l/2(p[j//](^))l/2 

1/ v®/l 


1/2 


< CB\f"{x)\'^ + C£B 


\F{x)\ ■ 


(4.13) 


Similarly, we have that 


T3,out<CBiV[r]ix))^/^ ( [ 

< Ce-^/‘^B^/^\f"{x)\^/‘^{V[f"]{x))^/^ 

< Ce-^B\f"{x)\^ + C£B^'£\^^. (4.14) 

\f [x)\ 


Combining (4.13)-(4.14), leads to the estimate 


2 \fix)T 3 ix)\ < Ce-^B\r{x)\^ + CeB^V[f]{x) (4.15) 


for some positive universal constant C > 0. This bound is consistent with (4.1). 
Estimate for the T 4 term. We decompose as 


r 4 = -8 P. 


'V-f 

Jr 

= 4/"(x) P.V. [ 
J\c 


(A„/(x) - f{x)){Aafix)f{5af'ix)f 


da 


da 


- 8 P.V. 

-sf 


((Aa/(x))2 + l)^a3 
{Aaf{x))^{f"{x) + ni[f"]{x,a)f 

H<V {{^af{x))‘^ + l)^ 

f (7^2[r](x,a))(A^/(x))2(^,f(x))(r(x)+7^l[r](x,«)) 

\a\<rj ((Aq/(x))2 + 1)^ Q;2 

{ni[f]{x,a)){Aaf{x)f{6af'{x)f 


da 


J\a\>ri ((Aq,/(x ))2 + 1 )^ 

= :^4, + 14,2,in + T4,out- 


-da 
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Using (2.24) we may estimate 
1 ^ 4 , 1 ,in I < C\f'ix)\ 


X / da + \f{x)\{V[f"]{x)Y/^ 


\a\<V 


\a\^^‘^da+ 'D[f'']{x) / \a\da 


'|a|<77 


< 


C\f W|DB|/ {x)\+eIBI +jB_^ 




< CeB\f'{x)y + C£^B^ 


\f''{x)\ 

By also appealing to (2.25) we obtain the bound 


|T4,2,in| < CB{V[f"]{x 


(If" 




(a:)l 


da 




\a 


11/2 


+ iV[f”\{x)f/‘^ / da 






On the other hand, for the outer term we obtain 

1 24.07.* I < CB^ 


[ 

J\a 


da 


<Ce-^B\f"{x)\ 


^|q;|>?7 I^I 

Combining (4.16)-(4.18), and using the Cauchy-Schwartz inequality we arrive at 

2 \fix)mx)\ < C£-^B\f{x)\^ + CsB^V[f]{x) 

for some positive universal constant U > 0. This bound is consistent with (4.1). 

Estimate for the T5 term. Recall that T5 is computed by taking the derivative of (2.3) as 


T^ix) = dr,v{x) = - 


= 2P.V. 


A 

dx 


P.V. 


1 


da 


(Aa/(x))2 + 1 a 


^af{x) daf{x) 


((A„/(x))2 + 1)2 «2 


da. 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


The main idea here is that in order to decompose into an inner and an outer term, we first need to subtract 
and add 


dd^,pv{x) — 2 . 


fix) 


,P.V. [ 

Jm 


daf{x) 


da = 2 


fix) 


a^ 


(/'(x))2 -p 1)2 

The contribution from this term is bounded by 

\n,pfx)\ < CB\Af{x)\ = CB\Hf"ix)\ 


(/'(x))2 + 1)2 


A/'(x). 


(4.20) 


where H is the Hilbert transform. The difference is now decomposed further as 

n - n,p, = P.V. [ Kf{x,a) (A„/(x) - /'(x))^^A/Mdce, 

Jm a 

where 

p- . 'I - n i^afix))^fix) + (Ac,/(x))2(/'(x))2 + /'(x)(A„/(x)) 3 + 2/'(x)A„/(x) - 1 

((A„/(x))2 + l)2(/'(x))2 + l)2 


Next we use that 


|iT/(x,a)| < 2, 


for any value of A„/(x) and f'{x). 
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Decomposing further, it is possible to find 

^5 — T^,pv = P-V- 


f , >(-a/"W/2 + K2l/"|(x.a)) (/»+Ki|/"|(i.a))^ 

— uajiwtw - S - 


J\a\>rj ((Aa/(x))2 + 1)2 


— • ^5,m “1“ 2^5,out* 

For the outer term we directly obtain 


\T5,out\ <CB^ 


i 


|/"(x)|. 

|a|>r7 1^1 ^ 


(4.21) 


We recall that rj = r]{x) = £B\f"{x)\ For the inner term, we appeal to (2.25) and (2.24) and obtain that 

\n,in\<\f"{x)\^ [ da + C\r{x)\{V[f']{x))^P [ \a\^Pda + CV[f"]{x) [ \a\da 

J\a\<v J\a\<v J\a\<v 


< C£B\f''{x)\ + C£^/^BP^{V[f'']{x))^P- 


1 


+ C£^B^V[f]{x)-. 


1 


(4.22) 


|/"(X)|V 2 1 /"(X)| 2 - 

Summarizing, (4.21) and (4.22) and using the Cauchy-Schwartz inequality, we obtain the desired bound 

\n{x) - n,p.{x)\ < c|ir(x)i + c£^b^'^^Q^ 

which combined with (4.20) yields the desired bound (4.2). □ 


Proof of Lemma 4.2. The proof closely follows that of Lemma 4.1, but uses (2.26)-(2.27) instead 
of (2.24)-(2.25). We let e E (0,1] to be determined, and as in (4.5) let 

V = V{x) = I 

\f ix)\ 

In this proof the constant C may change from line to line, and may depend on p, but not on e, U, or x. 
Estimate for the Ti,T 2 ,Ts, and T 4 terms. We first claim that 

|ri(x)| + |r2(x)| + |r3(x)| + |r4(x)| < ^\r{x)\^ + (4.23) 

for some constant C > 0, and all x E M. 

We verify this estimate by checking each term individually. For Ti, similarly to (4.6) and (4.7), by using 
(2.26)-(2.27) we have that 

\Ti,in{x)\ < Ce^/?’u(P+i)/P|/"(x)|(P-i)/P(DP[/"](x))i/?’ + 

^ d(p+1) /p 

SO that by Cauchy-Schwartz we obtain that Ti is bounded as in (4.23). 

For the T 2 term, as in (4.9), (4.10), and (4.11), but using (2.26)-(2.27), we obtain 

|T2,i,..(x)| < Ce-i/^>U(P+i)/?’|/"(x)|(^’-')/?’(P^’[/"](x))'/^’ + 

|T2,2,.n(x)| < cmf{x)\^ + 

|r2,„.t(x)| < ^|/"(x)|2 + ^!^|/"(x)|(^-^)/^(pp[/"](x))Vt’ 

SO that by Cauchy-Schwartz it follows that T 2 obeys the bound (4.23). 
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For T3, we proceed slightly differently from (4.6) and (4.7), but still use (2.26)-(2.27) and obtain 
ni[f]{x,a) A«/(x) 6af'ix) 


|7^3(x)| = 2 


P.V. 


\a 


|2-2/p ((A„/(x))2 + l)2 |a|2/p 


-da 


< CB{VP[f']{x)f/P 

< CB{VP[f"]{x)f/P 


(7^l[/'](x,a))P/(P-l)' 


(p-i)/p 


a 


(7^2[r](x,a)+a/"(x)/2)P/(P-l)' 


(p-i)/p 


'\a\<ri 


\a\ 




< 


CB{VP[f'']{x)f/P (|/"(x)|r?^/P + (pP[/"](x))i/V^P + 


which is consistent with (4.23). 

Lastly, for the T4 term, we use bounds similar to (4.16), (4.17), and (4.18), combined with (2.26)-(2.27), 
to deduce 

|T4,u„I < C|/"(x)| |^fB|/»(x)| 

\Ti2in\ < Ce^/PB(P+^)/P|/''(x)|(P-^)/P(pP[/''](x))^/P + 

’ ’ \f''[x)\PP 

\T,,out\ < ^\f'{x)\^ 

which concludes the proof of (4.23). 

In order to prove (4.3), we now use (4.23) in which we choose e depending on two cases: 


{B‘^/P{VP[f"][x)f/P < |/"(x)|2(P+i)/P, 


if 


then let 


e = 1, 

e-=|/"(x)|(i7PP[r](x))-V(p+i). 


\B‘^/P{VP[f"]{x)f/P > |/"(x)|2(P+i)/P, 

We thus obtain the desired estimate 

\T,ix)\ + \T2ix)\ + \T3ix)\+Ux)\ < C77max{|/"(x)|2,i72/(P+i)(pP[/"](^))2/(n+i)| 

< 


(4.24) 


CeB^VP[f"]{x) CB\f"{x)\‘^ 


+ 


x2/(p-l) 


|/"(x)|P-i 

where in the last inequality we have appealed to the e-Young inequality and the essential fact that 2 < p + 1. 
Estimate for the T5 term. The first step is to show that 


|r5(x)| < CB |'|A/'(x)| + i|/”(x)| 


(4.25) 


holds when p G (1, 2). For this purpose, similarly to (4.20), (4.21), and (4.22), by using (2.26)-(2.27) we 
arrive at 

\n,pv\ < CB\Af{x)\ 

\n,out\ < ^\f"{x)\. 
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Combining the above estimates proves (4.25). We now use (4.25) in which we choose e precisely as in 
(4.24) to obtain the desired estimate 

|T5 {i)| < CB (|A/'(i)| + |/"(x)| + B'/(r+i)(i,P[/»](^))i/(p+l)) 

< CB (|A/'(i)l + 4|/"(i)| + 

□ 


In the last inequality we have appealed to the e-Young inequality and the fact that 1 < p + 1. 


5. Proof of Theorem 1.2, part (i), Blowup Criterium for the Curvature 

We evaluate (2.7) at a point x = x{t) where the second derivative achieves its maximum, i.e. such that 

\nx,t)\ = \\r{t)u^ = M^it). 

At this point x we have 

dx\f''ix,t)\‘^ = 0 and Cf[f”]{x,t) > 0. 

By Rademacher’s theorem whose use is justified by the assumptions, we may then show that 

\\i^ = dt\nx,t)\^ 

< {dt + v{x,t)dx + Cf)\f''{x,t)\^ 

= 2f{x,t){T,+T2 + T3 + T4)-Vf[f]{x,t) (5.1) 

for almost every t G [0, T]. 

In order to complete the proof, we estimate the right hand side of (5.1). The dissipative term Vf[f"]{x) 
is non-negative, and is bounded from below explicitly as in (2.14), (3.1), and (3.3). In order to estimate the 
Ti,..., T4 terms on the right side of (2.6), we appeal to the upper bound (4.1) obtained in Lemma 4.1. We 
thus obtain 




d o 1 

—H- 

2(1 + 52) 


V[f"]{x) + \lb{Moo) < + C^eB^V[f"]{x) 


(5.2) 


for some universal constant Cq>1. We now choose the value of e in Lemma 4.1 as 

1 


£ = mm 


2Co52(1 + 52) 


,1 


in order to obtain 


where 


+ \lb{M^) < KbM, 


I 


Kb = CqB max{l, 4C7^5^(1 + 5^)2} 

is a constant that depends solely on B. To close the argument, we recall cf. (3.4) that 

,. Lb{M) 
lim ——^— = (X) 

M^-oo 

and that is continuous, which shows that there exists = M^{B, Cq) such that 


1 q 

-Lb{M) > KbM^ for any 


M>M^. 


Therefore, we obtain 

Moo(f) < max{Moo(0),M^} 
for all f G [0, T], which concludes the proof of the theorem. 
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6. Proof of Theorem 1.1, Local Existence 


In view of the available LP maximum principle for /, cf. (2.4), the proof consists of coupling the 
evolution of the maximal slope 

with that of the norm of /" 

M,{t) = \\nt)u.. 

Our goal is to obtain an a priori estimate of the type 

+ Mp) < polynomial(i?^ + M^) 

from which the existence of solutions on a time interval that only depends on i3(0)^ + Mp(0)^ follows by a 
standard approximation procedure. 

For the evolution of Mp{t) we split the proof in three cases: 

(i) p = oo 

(ii) p G [2, oo) 

(hi) pG (1,2). 

We however note that in all of these three cases, but the ID Sobolev embedding, as soon as f" G L^, we 
have /' G and we have the bound 

[/]c(p-i)/p < II/”IIlp 

for some positive constant C that may depend on p. In particular, it follows from (2.21) that the bound 

|7^l[/'](x,a)| = |A„/(x) - /'(x)| < \a\^-^'^/P\\f\\Lv (6.1) 

for all X, a G M, and all p G (1, oo]. In view of this estimate, we immediately notice that the uniqueness of 
solutions in follows directly from Theorem 1.4, whose proof is given below. 

6.1. Evolution of the maximal slope B{t). We first multiply equation (2.5) by f'{x,t) and arrive at 
the equation 

{dt + vd, + Cf)\ff + Vflf] = n (6.2) 

pointwise in (x, t), where we have denoted 

(7^l[/'](x, a))(Aa/(x, t)) Saf{x, t) 


Te{x,t) = 4/'(x,f) [ 
Jm 


rda. 


(Aa/(x,f))2 + 1 (h«/(x))2 + a2 

Let X = x(f) be a point such that \f'{x{t),t)\ = B{t). Using equation (6.2) and Rademacher’s theorem we 
find thaf 

^B‘^{t) = dt\f{x,t)\‘^ = Te(x,t) - (P/[/'](x,f) + £/|/'(x,f)p) (6.3) 

for almosf every t. Here we used fhaf af a poinf of global maximum dx\f'\'^ = 0. Nofe fhaf also thaf af x we 
have £/|/'P > 0. If remains fo bound Tq{x, t). For Ibis purpose we decompose Tg as follows 

f (B.i[f']{x,a)){Aafix,t)) 

'\a\<e 


Tq{x) = 2 
+ 2 


da 


(A„/(x, f))2 + 1 {6af{x)y + a2 
f (7^l[/](x,«))(Aa/(x,t)) 5a{f'{x,t)f 

\a\<s (A„/(x,f))2 + 1 ((5«/(x))2 + a2 

(7^l[/'](x, a))(A„/(x, t)) Saf{x, t) 


da 


((A„/(rc,i)P + l)2 

— 2^6,1,in( 3 :) “h TQ 2,in{x) T Tg ^^^((x) 


da 
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where e = e{x) > 0 is to be determined, and we ignore the t-dependence of all factors. Using (6.1) and the 
fact that Cf\f'{x)\‘^ > 0 we may estimate 

|r6,2,*n(x)|<e(P-l)/P||/"||Lp/:/|/'(x)|2, 


while the bound (2.20) yields 


Letting 


\T6,out{x)\ < 


£ = 


8B^ 


p/ip-i) 


\LP 


in the above three estimates, we arrive at 

\n{x)\ < \ (Py[/'](x) +£y|/'(x)|2) 


which combined with (6.3) gives 


A. 

dt 




3|| ^W||P/(P-1) 

Wlp 


for some positive constant C that may depend on p € (1, oo]. 


(6.4) 


6.2. Case (i), p = oo. We recall the evolution of Moo{t), cf. (5.2), in which we take 

e = eft) = min < , 1 > 

\2CoB{tY{i + B{tYy \ 

to arrive at the a priori estimate 

j^mI<CoBuib.^{2CoB\1 + B^),1]mI. (6.5) 

Combining (6.4) with (6.5) we obtain that 

< CB{tfM^{t) + CB{t){l + B{tffM^{tf 
dt 

<C{1 + B{tf+ M^{tyy (6.6) 

for some positive constant C. Integrating (6.6), we obtain that there exists T = r(||/g|| ^oo ^ \m ^oo ) > 0 
on which the solution may be shown to exist and have finite norm. 


6.3. Case (ii), 2 < p < oo. We consider the evolution of |/"|^ in divergence form, cf. (2.16), apply the 
upper bound given by Lemma 4.1 for the terms Ti,..., Ts on the right side of (2.16), and the lower bound 
given in Lemma 3.3, estimate (3.5), and the bound given by (2.14) for the dissipative term Py[/"] on the 
left side of (2.16), to arrive at 


I f"(x t)\P~‘^ 

{dt+cj) \nx,t)\p+dyv{x,t)\r{x,t)\p) + i^j-i_!i^p[/"](x,t) + 


|/"(x,f)|2P 


Cl(l + i?(t)2)||/'^ 


LP 




(6.7) 


pointwise in x and t. Choosing 


iCiB{ty{i + B{ty' 


■A 


£{t) = min 
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|/"(x,t)|2P 


we conclude from (6.7) that 

<Cii?(f)|/"(x,f)n77/"(x,f)|+Ci?(f)5(l + i?(f)2)|/"(x,f)ri 

where C = C'(C'i) > 0 is a constant. 

At this stage we integrate (6.8) for x G M. First we note that 

|/"(x)l^’-l/"(x-a)r 


p 

LP 


[ Cf[\fr]{x)dx = p.v. 

Jr 


rdadx = 0. 


( 6 . 8 ) 


(6.9) 


if{x) - fix - a))2 + a2 

This fact may be seen by changing variables x ^ x —a. We thus obtain an a priori estimate for the evolution 
of Mpit) = ||/"(f)||LP as 


—M (tY H_ ^2p{t) ^ _ 

dt ^ ’ Ciil + BitY)Mpit)P 


<CB{t)il + B{tffMp+i{t)P+^ 


( 6 . 10 ) 


by also using that the Hilbert transform is bounded on LP. Furthermore, since for p > 1 we have p +1 < 2p, 
we may interpolate 

Mp+iit) < 

which combined with (6.10) and the e-Young inequality, yields 

M2pitfP 


( 6 . 11 ) 


— M (t)P -\ _ 

dt ^ Ci(l + H(f)2)Mp(f)P 


< CBit)il + BitYfMp{ty-^M2pit) 


< 


M2pitfP 


+ CB{t)P/^P-^'>il + H(f)2)(3P+i)/(P-i)Mp(f)?’'/(P-^) (6.12) 


2C7i(l + H(f)2)Mp(t)P 
In conclusion, we obtain 

which combined with (6.4) gives 

d 


dt 


(H(f)2 + Mpitf) < C(1 + H(f)2 + Mp(f)2)5p/(P-1) 


(6.13) 


for some positive constant C. Integrating (6.13), we obtain that there exists T = T(||/q lli“) ll/ollip) > 9 
on which the solution may be shown to exist and have finite W‘^’P n W^’°° norm. 

6.4. Case (iii), 1 < p < 2. We proceed similarly to the case p > 2, but instead of applying to (2.16), 
we use (2.19). For the dissipative term on the left side of (2.16), we use the minimum between the lower 
bounds provided by (3.5) and (3.6) in Lemma 3.3 and (2.14). For the nonlinear terms on the right side of 
(2.19) we use the minimum between the upper bounds provided by Lemmas 4.1 and 4.2. The resulting 
pointwise in (x,f) inequality is 


idt + Cf)\f”\p + dMf''\n + 


1 


mm 


< C 2 eB^ min 


\f"?-P 

for some constant C 2 > 0. Choosing 




C2(1 + H2) 

C 2 B 




i/"p-p’ 


1/ 


n 1 2p 


C-2(1 + H2)||/"||^ 


LP 


e2/(p-i) 


\f'\P+^ +C2B\Hf''\ \f"\P 


(6.14) 


e(f) = min 


2C|H(f)2(l + H(f)2) 
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we conclude from (6.14) that 


(fl. + r,)l/T + 9A»I/T) + ^ 

< C(1 + B2)(p+7)/(2p-2)(-|j//|p+1 \ 



|jW|2p 

C2{i+B^)\\ni. 

(6.15) 


Upon integrating (6.15) for x G M, using the identity (6.9), the boundedess of H on L^, and the interpolation 
bound (6.11), we thus arrive at 


ity +__ 


< C(1 + 5(t)2)(P+7)/(2p-2)Aip^^(i)P+l 

< C(1 + fi(t)2)(P+^)/(2p-2)Mp(t)P-iM2p(f)^. 


(6.16) 


Similarly to (6.12)-(6.13), since p + 1 > 2, and by using the e-Young inequality, we conclude from (6.16) 
that 


d_ 

dt 


Mp{tf < C{1 + 5(t)2)(p'+9p-2)/(2(p-l)2)^^^^^^2^(3p-2)/(2p-2)^ 


Combining with (6.4) we finally arrive at 


^{B{tf + Mp{tf) < C(1 + B{tf + Mp(^)2)2p(P+l)/(p-l)^ 
dt 


(6.17) 


Integrating (6.13), we obtain that there exists T = T(||/g||p,cxj, ||/g ||lp) > 0 on which the solution may be 
shown to exist and have finite 11/2,p ^ pj/i.oo 


7. Proof of Theorem 1.3, Global Existence for Small Datum 

The proof follows closely the estimates in Section 6. The major difference is that assumption (1.4) and 
the maximum principle for /' established in [CG09, Section 5] show that 

||/'(f)||L- <i7< (7.1) 

Oj)s 

for all t > 0. Thus, we do not need to consider the evolution of B{t), as we have 

B{t) < T 

for f € [0, T), where T > 0 is the maximal existence time in lU^’^. For simplicity, we split the proof in 
three cases based on the value of p G (1, oo]: 

(i) p = oo 

(ii) p G [2, oo) 

(hi) pG (1,2). 


7.1. Case (i), p = oo. We use the estimate (5.2), but here we apply lower bound (3.1) instead of (3.12), 
and we set e = 1. We arrive at 

Mi, 


d ,-2 1 

-j\//2 _|- 

dt ^ 2{1 + Bi 


V[f]ix) + 


<CoBMi + CoB^V[f]{x) 


485(1 + 52) 

where x = x{t) is a point at which M^o = \f''ix, t)|. For 5 small enough, so that 


26*052(1+ 5^) < 1 


and 


10067o52(1 + 52) < 1 


dt 


^“^ 505(1 + 52)^-^^ 


< 0 . 


holds, we thus obtain 
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Integrating the above ODE we obtain that 


Moo(t) < 


Moo(O) 

1 I Mao (0) ^ 

^ lOOB '' 


for all t > 0, which proves (1.5). 


7.2. Case (ii), 2 < p < oo. We use the first line of estimate (6.8), but instead of using (3.5) to bound 
the dissipative term 'D[/"] from below, we appeal to (3.1). We arrive at 

< CiB\nx,t)\P\Hf"ix,t)\+CB^{l + B^)\rix,t)r\ 

Integrating the above over x € M, similarly to (6.10) we obtain 

1 


—M (t)P + 

dt ^ ^ ^2) + 4(1 + ^2) 


\f\x,t)\P-^V[f]{x,t)dx < CBil + ByMp+i{tr+^ 

(7.2) 


for some C > 0. If 77 is chosen small enough so that 

200CB'^{1 + B'^)^ < 1 

we thus obtain 




df'^P'-"' ' 2005(1 + ^2) 

Integrating the above in time and noting that 


+ 4(1 + ^V[f"]{x,t)dx <0. 


\f''{x,t)\P ‘^V[f']{x,t)dx = 


1 


(|/"(x)|P-2 + ir (x - a)|P-2)(/"(x) - fix - a)) 


-dadx 




4 

> - 
pZ 


{\f"{x)\Pl^-\f"{x-a)\Pl^f 


■dadx 






p- 


concludes the proof of (1.6). 


7.3. Case (iii), 1 < p < 2. We use estimate (6.15) in which we bound from below the dissipative terms 
from below by appealing to Lemma 3.1, and arrive at 

\f„\p+l 


{dt+Cf)\r\p+dMnn + 


+ 


2005672(1 + ^2 ) 672(1 + 52)||/"||P, 

< 67(1 + 5^)(p+'i')/(2p-2)(|j"|P+i _p \ \ f"\P). 

Integrating the above over x G M, and using that H is bounded on LP in this range of p, we thus arrive at 

A M (t)P H__^_ M2p{t) P _ j^2\{p+7)/{2p-2) (t)P~^^ 

dt ^ + 2005672(1 + 52) + 672(1 + 52)||/"||P, + ^ ^ ' 

Lastly, choosing 5 small enough so that 

400672675(1+ 52)(^P+®)/(2P"2) < 1, 

we arrive at 


(t)P H__ 

dt^ + 4005672(1 + 52 ) 


< 0 


which upon integration in time concludes the proof of (1.7) and thus of the theorem. 
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8. Proof of Theorem 1.2, part (ii), Blowup Criterium for Smooth Solutions 

We shall study the evolution of the ||/||/ffc(f) norms for k > 3. We show that they can be controlled 
by SUP[0 ,T] II/'WIIl- and sup[o,r] ll/"WII Then Theorem 1.2 concludes the proof. In fact, the 
norm of a solution with k > 3 can be controlled already by Tf^-norm as shown in [CCG+13a], Section 5.2. 
Therefore we may assume k = 3. We start by dealing with the evolution of ||/^^(t)||^2- We use inequality 
(7.2) with e small enough and p = 2, to obtain 


d ,, r/n, 2 I 

lli» + 


"l|2 
ijl/2 


2Co{l + B^) 


<C(i7)||/"||i3<C(i?)||/"|Uo 




and therefore 


II ^//||2 


I 2 exp (^C{B) ||/"(s)||Locds^ 


( 8 . 1 ) 


We use equation (2.1) to split 


f''{x)f”'{x)dx = h + h + h + h, 


where 


h= [ f"'{x) [ {f""{x)a-6o.f"'{x))( \ , 

7r 7r \{6af{x)y + a^ 

[ f'"ix) [ if'"ix)a - 5af''{x))dx f 
J M. J V 


I 2 — 3 


{6af{x)y + 


dadx, 

dadx, 


h = ‘i j f"'{x) j {f"{x)a-5af'{x))dl( ^ , 

7r 7r \{daj{x)y + a 

h = [ f''{x) I {fix)a - 6afix))dl ^ 

M IR. 


2 , dadx, 


{daf{x)y + a‘^ 


dadx. 


In Ii it is possible to decompose further and obtain 

h = [ rix)r'ix)p.v. ^ 

7r 




We bound Ii i as 


h,i = / \f"'{x)\^P.V. 


adaf{x)6af'ix) 


i{daf{x)y + a‘^y 


da 


= / irwpp.v. f !-( yy 

Jr jRa\{{Aaf{x)y 


fix) 


+ 1)2 ((/'( x ))2 + 1)2 


Aaf'ix)dadx 


+ 




2 f ix)Af{x) 

( X ))2 + 1)2 


dx + .fi,1,2- 


In one finds extra cancelation in such a way that splitting in the regions |a| < r and |a| > r and 
optimizing in r, it is possible to obtain as before 

h,i,i<ciB)\\ru^\\n\l2. 

For the Gagliardo-Nirenberg interpolation inequality IIsIIl- <C’ll9ll(7llslll7« yields 

/l,l,2 < ||/"||L^||/"'|li4 < C\\f'h4f'\\L4r'\\H^/2 


<C(i?)||/"||i.||/"'||i.+ 


11/ 


///||2 
H 


1/2 


32(1 + 52) 
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by using the e-Young inequality. This yields 


/i,i = < c{B){\\f"u<^ + iiriiioiir iii2 + 

Using (2.14), (2.13) together with (3.7) in I 12 we arrive at 

1 


f"' 11 ^ 

H 


1/2 


h,2 < — 


y //||2 

4(1 + .B2)"'^ "n 


32(1 + 52)' 


r"'ll3 


ll/"ll;i+ - 07/1 , r.2/3lir.ll II/^IL3- 


27(1 + 52)3||/"||^. 


Adding the last two estimates it is possible to obtain 


pf" 112 


nr iui /2 

h = h,i + h ,2 < c{B){\\r\\L^ + wfwhnrwh - 


n\h 


32(1 + 52) 27(1 + 52)3||/"||^, 

We are done with B. For I 2 we rewrite as 

^af”{x) - f"{x) Aafix) 


h =6 
=6 


/ fix) [ 

■/M -/M 

/ fix) [ 

7k 7|c 


a (Aa/(x))2 + 1)2 

7^l[/"'](x,a) A„/(x) 


l«l<ll/" 


+ 6 / /"'(x) 

t/ M. t 

•— B,in 4“ B,out- 

Inequality (2.24) allows us to get 




a (A„/(x))2 + 1)2 

T^l[f"]ix,o:) A„/(x) 

a (A„/(x))2 + 1)2 


Aaf'{x)dadx 
Aaf' {x)dadx 


Aaf {x)dadx 


l2,in<nf"\\L^ [ |/"'(x)| 5 [/](x)l /2 

7r 

<C\\f\\]!l [ \f'ix)\V[f]{x)^/^d: 

JR 


da 


UKiinizU 


dx 


<c{B)\\fun\n\h + 


II f"'ip 

32(1 + 52)’ 


and 


I 2 ,out B 125 [ |/"'(x)|5[/](x)V2 
7 k 


da 


l“l>IIT 


v/||-i |a|3/2 

r.oo I I 


dx 


<C 7 ( 5 )||/"||^/i [ |/"'(x)|P[/](x)V 2 da 

7k 

II f"'lP- 

<c'( 5 )iiriiic.|in| 2 ,+ 


32(1 + 52)' 

These last to inequalities give the appropriate bound for I 2 such that adding (8.2) we obtain 


r"'ii 2 


511/ iui/2 

/i+/2 < c(5)(iiriu^ + wfxnwrwh - 


n\h 


32(1 + 52) 27(l + 52)3||/"||^oo' 


( 8 . 2 ) 


(8.3) 
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It is possible to decompose further in I 3 as follows 


h = 


C 3 ,i [ f"'{x) [ (/"(x)a - 5„/'(x)) 
Jm. m. 




ii^afixW +a2)2 


dadx 


, f ttm ^ f tf"t ^ fii ^^{^af{x)Y{5af'{x)y 

+ C 3,2 f[x) if (x)a - Saf (x)) . — KT^dadx 


+ C 3,3 [ f"'{x) [ (/"(x)a-5„/'(x)) 

M. M. 

■=h,l + -^ 3,2 + -^3,3- 


{{Safix)Y + a'^Y 
daf{x)Saf"{x) 


i{Saf{x)Y + Q;2)2 


dadx 


The identity 


allows us to get 

^ 3,1 < 


f"{x)a — Saf'{x) = a^ f rf"'{x + (r — l)a)dr, 

Jo 

<C\\f"\\loo II I i\f"'ix)y + \f"'{x + (r - l)a)y)dxdadr 

Jo J\a\<B\\f"\\7L J^ 


+ CB^ 


/O JH>B\\f"\\iL 


{\f'"{x)y + \f'"{x + (r - l)a)y)dx^dr 


\a\ 


<cm\f''\\LAn\l^- 

An analogous approach for 13 2 gives 

42 <c{B)\\r\M\f'X 2 . 

For the 13^3 term, we decompose further: 

f'jx) - Aaf{x) A„/(x) 

'|o|<C 


-^3,3 —C3,3 


,3 / f'ix) I 
Jm. J Iq 


a (A„/(x))2 + l)2 


Aaf''{x)dadx 


+ C 3,3 I fix) I -4(/"4) - A fix)) 6a f"ix)dadx 

Jr J\a\><: ((Aa/(x ))2 + 1)2 


•—^3,3,in T lo,0,out' 

For the outer term, inequality 


yields 

Identity (8.4) together with 


\\f'\\h<c\\nLff\\L^ 


C*|l II „///!! , fill \\2 


h,3,out<j\\fY\ur\\L-<^\\n\h- 


A. 


J''ix) = I fix + (s - l)Q;)ds, 
Jo 


allow us to obtain for the inner term 

h,3,out < C3,3 I II I \fix)\\fix + ir-l)a)\\fix + is-l)a)\dxdadsdr 

Jo Jo J\a\<(;JR 

fW||3 


< 2 c 3 , 3 Ciiriii 3 . 


We take 


to find 


c = 


28 (l + i?2)3c3_3||///||^^ 

I 3 = h,l + h,2 + h,3 < CiB)\\f"\\L^ WfWl, + 




27(l + i?2)3||/n||^„,- 


(8.4) 


(8.5) 


( 8 . 6 ) 
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After comparing with Ii + I 2 in (8.3) we obtain 


5iini|v2 

/i+/ 2 +/3 < cmwru^ + iir iiioiinii^ - 


32(1 + .62) ■ 


(8.7) 


For the last term, we decompose using Leibniz rule to find 

h = 74^1 + /4_2 + 74,3 + Ia,a + 74,5 


where 


74,1 :=C4,i / f"'{x) [ {f'{x)a - 6af{x)) 
IR. »/ M 


Saf{x){Saf{x)f 
{{Safix)y + a^)^ 


dadx. 


T f f",, ^ f (ft ( ^ Si f( \S^ocf{x)f{bo,f{x)f 

74,2 :=C4,2 f [X) (/ {x)a - daf{x)) - Kzi-dadx, 


{{5af{x)f + a^Y 

Safix)6af"ix) 


h,3 :=C4,3 [ fix) [ if{x)a-6afix)) 

iz IR t/ M 

/4,4 :=C4,4 I /'"(X) /^(/'(X)« - &/(x)) ^ 


dadx, 


ii^afix))"^ + 
iSafix))^dafix)6af”ix 


dadx, 


h,5 :=C4,5 [ fix) [ if'ix)a - 5afix)) 

"ZIR </ M 

In 74 1 we bound as follows 
74.1 <C [ 


Safix)6afix) 


iiSafix))^ + a‘^)^ 


dadx. 


^ d«iiniL^ii/"iii8+cii/'iii.o [ 

7|a|<||/''||-i. J\c 

<ciB)\\ruff\\i., 


, iiriiL^iiriii.i^ 

I«I>II/"IIZ- 


where the last inequality is given using interpolation inequality ||/''||^8 <ll/'llL-ll/"lli^lirilL^ together 
with (8.5). The same approach allows us to conclude for 74^2 that 

74,2 <C(77)||/"||i^||/'"||i,. 

Using identity ( 8 . 6 ) in 74^3 we arrive at 

"1 


74,3 <cii/"iii.oiinii. 


'0 7|a|<||/"||-i, 

fW||2 


dads + C\\f\\f\\f\\l^ n 

Jo J\c 


da 


0 7|a|>||/"||-i, |a| 


ds 


<CiB)\\f\Lff\\f 

The same procedure yields 

hA<ciB)\\n\Lff\\f 

Finally, in 74^5 we use next splitting 

74,5 <C||/"||l^ [ \fix)\( [ 

7r ^J\c 


da 


i-i 

IliOO 


1/2 


O N 1/2 

ix)\^da] dx 
da 




+ C\\f\f [ 

J\c 


\a\>B\\f" 


H<B\\f" 

ilfif"^ + \fix - a)\'^)dx 


\a\ 


<CiB)\\f"\\^/l / \fix)\iVif)ix))^/^dx + CiB)\\f'\\Lff\\l2 


<C(4B)||r||i^||rili, + i|^l^. 
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Above estimates allow us to conclude that 

5 

/4 = ^/4,fc<C(i?)||/ 


> 11 / 


W ||2 I ll/^”llijl /2 


+ 


k=l 


Adding to (8.7) we obtain that 


^ " fWll2 , ll/^^^lliTi/z ^ .f"ll2 


dt 


iir 


lL2 


+ 


8(1 + 52) 


' 32(1 + 52)' 

<c(5)(ii/"iii. + ii/"iii^)iiriii.. 


The bound for ||/”||l 2 in (8.1), the control of 


and integration in time yield the desired result. 


9. Proof of Theorem 1.4, Uniqueness 

We consider two Muskat solutions fi and /2 satisfying the hypothesis of the theorem with the same 
initial data fo{x). From (2.1) and a small computation we obtain the equation for the difference g = fi — f 2 , 

dtg + vidxg + 5/^ [g] = T 7 (9.1) 

where vi{x, t) is as defined in (2.3) in terms of /i(x, t), Cf^ is defined as in (2.8), and 


T-r{x,t) = P.V. 


Sag{x,t) {f2{xp) - Aaf2{x,t)) {Aafl{x,t) + A«/2(x,f)) 


da. 




(1 + (A„/i(x,f))2)(l + (A„/2(a:,t))2) 

Let 5 = supjg[o,T]y=i ,2 ll/jiloo- By assumption, has a uniform modulus of continuity p on [0,5], and 
thus by (2.21) we may find an e = e(5, p) > 0 such that 

\f 2 {x,t) - Aaf 2 {x,t)\ < p{\a\) < p{e) < ^ (9.2) 

for all jaj < e, and all (x, f) G M x [0, Tj. We fix this value of e throughout the rest of the proof. Denote 

{f 2 ix,t) - Aaf 2 {x,t)) {Aafl{x,t) + A„/ 2 (x,f)) 


Ki^2ix,a,t) = 

It follows from (9.2) that 


(1 + (A„/i(x,f))2)(l + (A„/2(x,f))2) 


|A'i,2(x,a,f)| < ^ 

for all I a I < e, while the Lipschitz assumption on /i and /2 directly yields 

\Ki^ 2 {x,a,t)\ < 25 

for all I a\ > e, uniformly in x and t. 

Upon multiplying (9.1) by g{x, t), and recalling the definition (2.9), we obtain 

{dt + vidx + /i/Jjpp + Pfiig] 


(9.3) 

(9.4) 


= P.V. 


5a{g^{x)) 




7|o|<£ 

+ 2g{x, t)P.V. 


Kipix, a)da + P.V. 




((5„p(x))2 

-^-ATi 2 (x, a)da 

1/-^ ’ 




^ag{x) , . 

- 2 —-^ 1 , 2 ( 3 :, a)da 


J\a\>£ ^ 

= ■ T 7 ,l,in + 57^2,in + Tj^out- 

First, we notice that in view of (9.3) we have 

1 57,2,in I < 

while in view of (9.4) and the Cauchy-Schwartz inequality, we have 

\T,,out\ < 45|<7|(5yj5])V2e-i/2 < ^Vf,[g]+2B^e-^\9\". 


(9.5) 
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The above two inequalities combined with (9.5) yield that 

{dt + vid^ + Cf,)\g\'^ < 2B^e-^\g\^ + (9.6) 

To conclude, we note that the decay assumptions of the theorem guarantee that there exists a point, denoted 
by X = x{t), where \g{x,t)\ = || 5 r(f)|| At this point of global maximum we have that dx\g\'^ = 0, and 
> 0. Moreover, 

\T-J^l^in{x{t),t)\ < ^(£/Jpp)(x(f),f) 

and thus from (9.6) we obtain that 

{dt\g\‘^){x{t),t) <2B‘^e-^\\g{t)\\\o.. (9.7) 

The pointwise differentiability assumptions further warrant the use of the classical Rademacher theorem 
(see Appendix) which implies that 

^Il9(*)lli~ =(ail9l")(3:(*),*) (9.8) 

for almost every t, where x = x{t) is as above. From (9.7), (9.8), and the Grdnwall inequality it follows that 

||9(i)l|L°° < lb(0)||L°° exp {2B‘^e-^t) 
which concludes the proof of the theorem since ^((0) = /i(0) — / 2 ( 0 ) = 0. 


Appendix A. Rademacher Theorem 


Let us recall the classical Rademacher theorem for the convenience of the reader. Suppose /(x, f): M x 
M+ —)> M is a Lipschitz in time function uniformly in x, i.e. 3L > 0 such that for all t, s, x we have 

- f{x,s)\ < L\t - s\. 

Suppose that at any time t there is a point x{t) such that f{x{t),t) = M{t) = sup^, /(x, t). Then M{t) is a 
Lipschitz function with the same constant L. Indeed, let t,s £ M+ and wlog M{t) > M{s). Then 

M{t) - M{s) = f{x{t),t) - f{x{t),s) + fix{t),s) - f{x{s),s) < f{x{t),t) - f{x{t),s) < L\t - s\. 

V 

<0 


Then by the classical Rademacher Theorem in ID, M is absolutely continuous on any finite interval, i.e. 

M{t) — M{s) = f m{T)dT, 

J S 

where ||m||oo < L, and hence M' = m a.e. To show that M'{t) = dtf{x{t),t) at the points where M' 
exists we need extra assumption: for every x, /(x, ■) is differentiable everywhere in t. Then 


M'{t) 


lim 
h —^ 0 -|- 


f{x{t + h),t + h) - f{x{t),t + h) + /(x(f),f + h) - f{x{t),t) 

h 


> 


lim 
h —^ 0 -|- 


/(x(f),f + h) - f{x{t),t) 
h 


dtf{x{t),t). 


Taking h < 0 proves the opposite inequality. 
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